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Abstract 

The evolution equation for high twist operators is solved in the 
Double Logarithmic Approximation for cylinder-type diagrams which 
dominate in the limit of large number of colors. The asymptotic be- 
haviour of the parton correlation function is shown to be determined 
by the spectrum of the Ising model for a one-dimensional chain of 
spin-1/2 magnets. 



1 



1. The strong growth of the deep- inelastic structure function with de- 
creasing Bjorken variable xb leads to the necessity for taking into account 
the contribution of high twist operators since when xb —>■ the power sup- 
pression l/Q^ can be compensated by the large logarithmic factors ~ IuI/xb- 
Indeed, the moments of the structure function 

= dxBxi-'xBG{xB, Q'), (1) 
Jo 

where is a square of virtual photon momentum, can be expressed through 
a Wilson Operator Product Expansion in the form (see[0]) 

M{j,Q') = E ^CMj,Q'){h I 0(2^) I h). (2) 



N>0 



Here C2Ar(j, Q^) is the coefficient function and {h \ O^"^^^ \ h) denotes the 
expectation value of the local operator in a hadron state. Usually one can 
neglect the contribution of high twists (i.e. the terms with > 1 in (^) 
as they are suppressed by the powers of the large momentum square Q^. 
However this is not the case for ^ 0. As is known from renormalization 
group arguments the coefficient function is expressed through anomalous 
dimension 72jv(j — 1) 

C2iv(j,Q')~e^-(^-^)''^«'/'^', (3) 

where /i^ is an infrared cut-off. The small-xe behavior of the structure func- 
tion is controlled by the singularity of 72Ar at j = 1 . For the twist 2 case 
the anomalous dimension is given by the kernel of Gribov-Lipatov-Altarelli- 
Parisi (GLAP) equation in the gluon-gluon channel, and near j = 1 it has 
the form: 

720-l,Q) = ^^(— ^ (4) 

[Nc is the number of colors). The quark contributions are negligible here 
since they are not singular for j ^ 1. Supposing the high twist anomalous 
dimension '-)2n to have the same singularity at j = 1 one can estimate the 
contribution of the twist 2A^ operator to the structure function for a;^ ^ 

as 

/ 9 \ 



xbG2n{xb,Q^) ^ exp A2Aryln(52//^^lnl/xij, (5) 

where the coefficient A2Ar is determined by the residue of 72Ar at j = 1 .This 
relation suggests existence of a kinematic region of the xbiQ'^ variables where 
the power suppression is compensated by the logarithmic factors in the ex- 
ponent. The line where the twist 2 and twist 4 contributions are of the same 
order is the boundary of this region, and in the interior of it the twist 4 
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becomes leading. There is an infinite set of similar domains where the twist 
6 operators turn out to be dominant ones, then twist 8 and so on. The exact 
location and boundaries of these domains are determined by the coefficients 
A27V in Eq.(^. Inside these domains GLAP evolution equation becomes invalid 
and has to be replaced by the appropriate equation for high twists. 

In paper the specific contribution to the anomalous dimension '~f2N 
coming from the exchange of 2N non-interacting ladders (Pomerons) was 
found to be 

72jvH =iV72(f), ^ = j-l (6) 

or 




xbG2n{xb,Q'^) ^ [^j expN]j^^^\nQ'^/fi'^\nl/xB (7) 

that is the product of one-Pomeron asymptotics. 

The calculations for twist 4 case performed in Refs.p, |^ with the ac- 
count of Pomerons interaction in Double Logarithmic Approximation (DLA) 
modified this result: 

74M = 272(^)(l + 5) = ^^(l + 5) (8) 

Z ITU! 

with 6 ^ {N^ — ly ~ 10^^. The corrections due to other states (color lad- 
ders) were estimated to be approximately of the same order. The 6 value 
which indicates the deviation from Eq.@ can be treated as a measure of 
the Reggeon interaction strenght. The problem how the Reggeon interaction 
modifies result (|^) will be concidered below. 

2. This paper will be based on the direct solution of the evolution equation 
for high twist operators in DLA. 

The evolution equation for the high twist operators was derived in paper 
[^]. It generalizes the GLAP equation for quasipartonic operators which form 
the closed set of high twist operators allowing for the interpretation in terms 
of the parton model. Decomposing the momenta of virtual gluons ki into 
Sudakov variables 

ki = aiq' + Ap' + k±i, {p'f = {q'f = 0, k±ip' = k±,q' = 0, 

where p' is a hadron momentum, q = q' — Xbp' is a photon momentum, 
q^ = —Q^, matrix elements of the twist quasipartonic operators between 
parton states is expressed only through the jSi variables and has a general 
form 

pS:;;;:r..(A, • • • ,/?^) = r;:v;.:;-/?r ■ ■ -/^r, (9) 

where are integer numbers and tensor F specifies the color and Lorentz 
structure of the operator. The variable /3j has a meaning of the fraction of 
the total hadron momentum carried by a parton. 
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In axial gauge g^A^ = only the ladder-type diagrams contribute to the 
leading logarithms. For the twist case they comprise the local operator 
vertex of the form (|^) and gluons in t -channel interacting through all 
possible s-channel gluons rungs. Very important property of Leading Loga- 
rithmic Approximation (LLA) is strict ordering of transverse momenta along 
the ladder graph, namely, the transverse virtualities decrease from the top 
of the diagram to the bottom. The range of integration in each ladder cell 
is such that transverse momentum above the rung is greater than the mo- 
mentum below it. In other words the above momentum plays the role of 
ultraviolet cut-off for the below cell. Thus the value, being the great- 
est momentum in the upper loop that is attached to the local operator, is 
the ultraviolet cut-off for the whole graph. The evolution equation for the 
quasipartonic operators is derived from the diagrams with pure longitudinal 
external gluons' momenta Pi ~ f3ip'. Their sum determines the parton corre- 
lation function p{Q^, ■ ■ ■ , Pn) the evolution equation for which is obtained 
by taking the derivative of the diagrams with respect to InQ^. It has the 
form of the A^-particle one dimensional equation with pairwise interaction in 
which the variable ^{Q^) = /^f ^^^^ is an effective "time" ias{k^) is the 
running coupling constant) 

-p(e, A, ■ ■ ■ , = 5: / dp'.dp'^ 5{P, + Pj - A' - X (10) 

x$(A,/3, I A',/5;)p(e,/3i, ■ ■ ■ , A', ■ ■ ■ ■ ■ ■ 

(to be more compact the color and Lorentz indexes are suppressed here). The 
kernel \ (3'^, Pj) is determined by the logarithmic part of the one-loop 

integral over k± and a. The initial condition for the function p(^) is given 
by the expression (H) that corresponds to the absence of loops. The equation 
( p!OD collects all the powers of as In Q"^/ p^ where the infrared cut-off /i^ -C 
is the characteristic hadron scale. 

The longitudinal momenta /3j are not in general ordered in LLA. But in 
DLA they also have to be ordered to provide a large logarithm for each ladder 
cell. The Pi variables increase from the values ~ at the local operator 
vertex to the order of unity ones in the lower part of a diagram. In such a 
kinematic the logarithmic divergencies that occur in every loop when P'^ —>■ 
is cut from below by a longitudinal momentum in the upper cell. Thus DLA 
implies the loop integration in the evolution equation ([ToD to be limited by 
the condition 

(11) 

which means that the momenta below the s-channel rung (P) and above 
it [P') are of different order of magnitude. The most singular contribution 
comes in Eq.(|lOD from the region where both Pi and P'j tend to zero. The 
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momentum conservation allows it only if 

A + (12) 

that is with logarithmic accuracy I3i ~ Hereafter it is convenient to 

assume the momenta directed upward to be positive, directed downward to 
be negative. 

In the logarithmic domain the kernel of the evolution equation can be 
easily obtained in helicity basis [^] by keeping the most singular in (3' terms 
in gluon-gluon channel and integrating it over and a 

<.fAi(A,/?2 I = 25,,.5x,.ir''^Hf^'^'mPi-P2) (is) 

P1P2 

Here /x, z/, A, a are the two dimensional transverse indexes , f°-''^'3 are the 
structure constants of the SU{Nc) group. The momenta /?i,/?2, are equal in 
DLA . They are positive but have the opposite directions, one of them is 
incoming from below the other is outgoing. The low-scale momenta f3[ and 

are not supposed to be equal in DLA since the momentum transfer from 
below Pi — P2 is small only compared to the large momenta /3i and P2 but is 
of the same order as the low-scale ones. However P[ and P'2 like Pi and P2 
have the opposite directions, so the incoming from below momentum turns 
out into outgoing and vice versa. This property of the longitudinal momenta 
directions can be formulated as an "arrow rule" 0. One can assign an arrow 
to each t-channel gluon line to indicate the momentum direction. The inter- 
action occurs only between the lines with the opposite arrows while for the 
particular line the arrow direction does not change through all interactions. 
As a consequence all t-channel lines are divided into those incoming the local 
operator vertex and outgoing it. The violation of the "arrow rule" results in 
the loss of the longitudinal logarithm in the cell integration producing only 
small factor ~ as- One should note here that the "arrow rule" will be ef- 
fectively absorbed into restrictions imposed in the Nc ^ 00 limit where only 
the "neighbouring" gluons are allowed to interact. 

It should be noted that the restriction only to the ladder type diagrams 
is valid in DLA provided the colorless in t-channel states are considered. In- 
deed, for these states only hard s-channel gluons carrying the longitudinal 
momentum fraction P^^^ = Pi — P[ ~ Pi are emitted. The logarithmic con- 
tributions from the soft non-ladder gluons emission cancel out when all the 
diagrams are summed up, since the soft particle probes the total color charge 
of the system which is zero for white operators 

The kernel ( [TB| ) gives rise to DLA evolution equation 

d 
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X 



k<j 

fPk d(3'^ ff^^ dPj au...,bk,...,bj,...,aM(f) f)' f)' /Q\ 

lo Jo /'mv..,m;^iv..,<Tiv yPi^ ■■■iPki---iPji---iPN)- 

Here all (3i > 0, the momentum direction and the "arrow rule" are specified 
by the indexes a and coefficient S^^, ^„ . . The integration range in Eq. ([l^) is de- 
termined with logarithmic accuracy by condition ( |1T]). The initial condition 
is given by Eq.(^ 

For twist-2 colorless operators the solution to Eq.(|l^) has the form 

where function / obeys conventional GLAP evolution equation after substi- 
tution P = I (the structure function is considered to be dependent on the 
variable -77- = ^). 

2pq 13 ' 

The color structure of the evolution equation is essentially simplified in 
the Nc ^ 00 limit. If in this case the color index is written as the pair of 
quark and antiquark ones a = (i,j) and gluons are drawn by double lines the 
interaction occurs only between the gluons containing the common "quark" 
line. The closed cycle arising when such t-channel gluon lines are connected 
produces a factor ^Nc while joining other lines gives an order of unity factor. 
There is a significant difference here with the twist-2 (Pomeron) solution for 
which the factor Nc is obtained after adding together two possible joining of 
two t-channel gluon lines. The gluons can be enumerated in such a way that 
the gluon with number n is allowed to interact only with the gluons n — 1 and 
n + 1. All the gluons can be represented to be lying on the cylinder in the 
color space and each of them interacts only with its neigbours. It is natural 
to assign to the neighbouring gluons opposite signs of jS to ensure interaction 
between them. After this the color structure of the function p turns into 
itself under the action of the kernel of the evolution equation which takes the 
form 

d 

^Pmv-,/^]v(A) • • • ^Pn) = 

= II'^Mn,Mn+l'^<K+i f^nSiPn - Pn+l)y< (15) 
n 

^0 JO Pn+1 

Here is redefined (for large case) as 
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3. There is a way to transform evolution equations (0) or ( |T5| ) to the 
Schrodinger-type differential equation. The Laplace transformation 

PCX) 

p{E) = / rfee-^«p(0 (17) 



brings the evolution equation to the stationary form, which in the logarithmic 
coordinates Xj = In reads 

^Pri,...,riv(a;i) • • • ^a^Tv) = (18) 
/■oo roo 

= S{xi Xk) / dx^ / (ix^ j4jfc pr-^^,,,^j.jy (xi, . . . , Xj, . . . , Xjij, . . . , ajjv)- 

The arrow,Lorentz and color indexes are accumulated here in the index 
r = {a, /i, a} and the matrix Aik acting on the indexes at i, k locations is 
introduced 

A,, = 25.^,.,5,^,,,5,, z/'^''<'^^/^'<''^'=. (19) 
Seeking the solution of Eq. (|1^) in the form 

d d 

Pri,...,rjv(2^1) • • • ) ^N) = ■ ■ ■ V'ri,.,rjv(2^l5 • • • ) ^A^) (20) 

the function ip obeys the equation 

N g 

E'[[-g-'^r,,...,rj,{Xi,...,XN) = (21) 

i=l 

d 

= -'^S{Xi-Xk) n -K— ^ik<^ru-,rNi^ly ■ ■ y^N) 

{ip is assumed to be vanishing at the upper limit). The coordinate space is 
separated into A^! sectors which are obtained from the fundamental one 

xi < X2 < ■ ■ ■ < xn-1 < xn (22) 

by the permutation of any two variables. The function ip satisfies inside these 
sectors the free equation 

N Q 

^V''-i,...,r-jv(2^i' • • • ^a^Tv) = (23) 

and has the form 

<^,,„..,,^(xi, . . . , x^) = e*^^^i+-+^-^- (24) 
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or linear combination of such terms with some ki equal to zero in each of 
them. At the boundary separating the two sectors the matching condition 
has to be fulfilled. Near, say, xi = X2 hyper plane it reads 

^ d d d ^9 

Denoting the solution in the sectors Xi > X2 and Xi < X2 as and ip"^^ 
respectively, so that 

if = e{xi - X2)ipll^^^^ (Xi, ...,Xn) + 0{x2 - xi)ip'^^l^^^^^^{xi, ...,Xn) (26) 

the functions ip are connected through the relations 

2) E [(^^- - afc^"T+ (A^" L,-.. = (27) 
= -A,2{r' + ^'')l,=^,. 

Conditions ( pTl) allow to establish the connection between the evolution 
equation and one- dimensional quantum mechanics of particles described 
by the Schrodinger equation 

E 92 

-^Y^g-^ <^r-i,...,r^ - ^i^i - ^k) Aik V?ri,...,r^ = S^r,,...,rj, ■ (28) 
j j i<k 



Indeed the function (p clearly has form (|2^) inside any sector of type (g^) and 
satisfies conditions (^) at its boundary. Thus Eqs.(^) and ( pHj ) have the 
same wavefunction although there is not straightforward connection between 
the energies E and S. For the large Nc case Eq.(PBD can be rewritten in a 
chain-like form 

E 

(xi,...,xn) = Sip^,^,„^f,^{xi, . . . ,xn), (29) 

^ n 

1 92 92 

The correlation function p built upon function ip through (^) represents 
the sum of the terms with various number of 5-functions. Note that the term 
without 5-functions vanishes after some ki are put to zero. Its absence is 
evident also from the structure of the equation (pT|). 

4. Instead of the solution of evolution equation (|T^ or (^) the approach 
based on the direct summation of the ladder-type diagrams in the ^ oo 
limit will be adopted here. To begin with the formal solution to the evolution 
equation (|1^) can be written as 

00 1 

p(0 = e*«Po = E-e"<^'>, (30) 

n=0 ^■ 
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where $ is the kernel of the evolution equation. The term acting on the 
initial state po generates the ladder-type diagrams with n loops in /3-space 
(with n integrals over 13). To restore the proper expression (|30|) each /3-loop 
has to be multiplied by variable ^ and the whole diagram has to be divided 
by n\ . Introducing function 

n=0 ^ 



the solution (|30|) is given by the formula 

dE 



where the integration goes aroung = point in the counter-clock-wise 
direction. Function W{E) obeys Bethe-Solpiter equation 

W = + ^^W (31) 

describing the system of t-channel particles with pairwise interaction be- 
tween particles n and n + 1 

^,,u,xA(3i.(32 I = S,,Jx,.PiS{Pi-P2) (32) 

P1P2 

Equation ( pT] ) gives rise to the set of ladder-type planar diagrams in /3-space. 
To sum it up they can be transformed to the form where the interaction 
between t-channel gluons is treated as an exchange of s -channel particles. 
Kernel (0) results effectively from the exchange of the particle whose emis- 
sion is given by the expression (Fig.l): 

T,. = -^0{(3-(3')^^al{(3), (33) 



where aj^i^{P) {axa{j3)) is a creation (annihilation) operator 
axaiP) , aluiP') = SxfiSavS{P - (3'). 



Indeed, amplitude (^), being squared, reproduces evidently (|32D . 

In a ladder-type diagram any t-channel line emits and absorbs arbitrary 
number of s-channel particles. For the particular t-channel line with number 
n there is associated an amplitude containing creation operators cL}^^^y{f3) and 
annihilation ones an-i,\a{l3)- For instance, the amplitude shown in Fig. 2 has 
a form 
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^ V 2; -py ypi ^ \/P2 ^ X 

Pi P2 Ps 



/3 V '-4" 

P4 ^ y 

where the integration over the longitudinal momenta is submitted to the 
restriction 

x> pi> p2> f33> P4>y- (34) 

Multiplying expressions of this type and applying Wick theorem to the prod- 
uct of the creation and annihilation operators one can reproduce all the 
ladder diagrams for the equation (|3l|). Note that only the planar graphs 
will appear since the diagrams with the crosses between s-channel particles 
destroy (3 ordering (|3^) . Assembling for a given t-channel line the terms with 
all possible emissions and absorbtions one gets the expression 



X ^ (an,^,^,^{x) + a|,_i^^^^ (x)) ^ (vmim^IA) + an-i,Mi/.2(A)) x (35) 

■ ■ ■ X ^ (an,M,M'(/3fc) + an-l,Mfe/.'(A)) • 

Here x and y are the momenta of lower and upper points of the t-channel line 
(incoming and outgoing momenta, see Fig.2), /i and 11' are their polarization 
indexes. Introducing 2x2 matrix An{l3) 

An,^,A(^) = a„.,.(/3) + 4-i,M.(/5) 
Eq.(|35|) can be rewritten as 



VWEy ' Jy '^Jy Jy ' VWl VWk 



X 



Xb 



A„(a;)Pexp / d(3^^)\ . (36) 



En, 



Here symbol P means ordering of v4„(/5) matrixes according which the matrix 
with larger argument stands to the left from the matrix with smaller one. 
An equivalent form of the latter expression is 

-/■^<--)-^^|(--{//''^}) (37) 
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Now the problem reduces formally to the evaluation of the correlator 



Mi,A.i,...,Miv,MW(a;i,2/l, • ■ ■,XN,yN) 











(38) 



where |0) is the vacuum of operators ai, a2, . . . , Qn- Indeed, only Wick pair- 
ing terms survive after vacuum averaging, therefore correlator (|3^) includes 
all the ladder diagrams for Bethe-Solpiter equation (pT]) , whose solution is 
expressed through it as 



W, 



(Xi, . . . , Xn) 



■ dvN 



/II, IM, 



For the cylinder topology case one has to replace in (|38| ) aj — > a^r in the 
first term and oat — ajy in the last one. It is needed to ensure Wick pair- 
ing between t-channel lines 1 and N to reproduce their interaction on the 
cylinder. 

If it were not for the matrix structure of T^"'^ the correlator would be 
easily calculated using the relation 



exp{ / d(3U{f3) anim exp{ / d(3V{(3) at (/5)} 



(39) 



= exp{ J d(3V{(3)al{(3)} exp{ J d(3U{(3) a„(/?)} exp{ J df3 U{f3)V{[3)} 

which is valid for arbitrary functions U{f3) and V{f3). Indeed, in that case 
the correlator would be of the form 



where 



,XN,yN) 



N xi xn d 

E2 



d 



yi 



Hn dxi dx 



■ip{xi,yi, . . .,XN,yN), 



N 



exp 



X exp 



yi 

^2 dp 
y2 



(p{xi,yi . . .,XN,yN) 
ai{fJ)> exp 



=02 (/3) \ exp 



yi 



dp 

7W 



X 



(40) 



a[{p) 







The evaluation of this expression can proceed as follows. To begin with one 
can get rid off the fourth exponent in (^) moving it to the left through third 
and second ones, commuting it through relation ( |5D| ) with the first one and 
taking into account that (0|a|(/5) = 0. Then the same can be done with the 
other exponents moving the terms with creation operators to the left and 
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with annihilation ones to the right. After a sequence of such manipulations 
the final result will be 



ip{xi,yi,...,XN,yN)= (41) 
_ I /-^la dp ^ r^^-i dp ^ ^ r^N~i,N dp ^ r^N,! dp 1 

\Jyi2 EP Jy23 EP JxN-l,N EP Jxn,1 EP J 

where the integrations limits Xik = min{xi,Xk} and yik = max{yi,yk} arise 
after substitution in relation (||) Ui{p) = Vi{p) = ^= e{xi - P)e{p - yi). 



Differentiating these terms produces pre-exponent factors 6{xi^i — Xi). Their 
presence in the correlation function has the origin in the evolution equation 



structure (see Eqs.(^) and (26)). In DLA framework all Xi in the exponents 
have to be taken to be the same and all y^ c:^ xb- Then the leading behaviour 
of the correlator will be 

(O T'-^\x,xb)---T^''\x,xb) o\ ~exp{Arlln— }. (42) 

\ ' E xb 

5. The matrix structure of (^) prevents to apply relation (|39|) immediatly. 
To avoid this difficulty P-exponent in (^) can be represented through the 
functional integral. This expression is derived in two steps. At first the P 
interval is divided into large number N of small intervals A/?: 

P expj r -^A(/5)| ^ (43) 



/Ep 



V VEPn J \ ^EPn^i ) \ y/EPo 

pN = X Po = y, 

the accuracy being the better the smaller value A/3 is taken. At the second 
step the following expression for the unit 2x2 matrix is inserted between 
the terms in the product 



c J dz*dz 5{z* z — 1) ZxZ* = 5xa- 



Here c is a normalization constant and the integration is carried out over the 
components of two dimensional complex vector z, 



dz*dz = dz2dz2 dz\dzi 
5{z*z - 1) = 6{zlz2 + zlzi - 1) 



It gives 



( r dp ] /•^+^ 

^expM^ -=A{p)\ J l[ dzldzk6{zlzk ~ I) zn+1 
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Denoting = z{l3k) and taking into account that for = 1 

the product in (^) takes the form 

fc=o V V^A^fc / 

~ exp [ {z*z + z*-^£=Az)d(3, 
JxB V Ep 

where dot means the derivative over (3. As a result the amphtude is given 
by the functional integral over the complex vector field z{l3): 



PilXx,y) = c' I Dz*^DzJ{z*^Zn - 1) Zr,,^{x)z*,{y)x 



X exp <^ I dp 



(45) 



Here d is the normalization constant and 

Dz*^DZn6{zl^Zn - 1) 



: xYdzmdznim^m^niP) - 1)- (46) 

Instead of an explicit (5-function in measure (^6]) one can rewrite the integral 
as 

Ptl^'{x,y)=c" J Dar^J Dz:Dz^x 
d(3 



X exp 



(47) 



where auxiliary field (Tn(/5) is introduced. The functional representation (^) 
or (|47|) has to be substituted in correlator (|38|) . Averaging of exponential 
factors over vacuum state can again be done by making use of formula (p9|) 
taking 
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It results into expression 



(48) 

I PXr) PXn n -I- 1 1 

iVn) F{z*, z^ exp <^ 



y~z„ 



where pre-exponential factor F{z*, z) arises after differentiating the exponent 
over Xi variables: 

F{z*, z) = ]X[xizti^i)^ii^i) + hi_i^i9{xi-i - Xi) + hi^i+i9{xi+i - Xi)]. (49) 

i 

Prime here denotes that after multiplication one has to replace 

{hik9{xi - Xk)f hikXi5{xi - Xk). 

Within DLA accuracy all Ui in the exponent have to be taken to be of the 
order of xb while all ^ a; ~ 1. After that integral (^8]) can be treated as a 
functional integral for a one- dimensional quantum mechanical chain. 

Let X and Cn,\ be the creation and annihilation operators of the particles 
labelled by the index A in the chain site n and let the dynamics of the system 
is described by the hamiltonian 

n 

which commutes with the operator of the number of particles in site n 

2 

A=l 

Writing (5-functions in the measure of Eq. (^8]) through additional integrals 
over functions cr(/3) the integral over z*,z variables can be considered as a 
functional integral in the holomorphic representation written for the matrix 
element 

M = (50) 



0. 



Y[cm,,^Fic\c) exp I r rf/3Eza„(/3)(Ar„ - 1) + 1^H]]y[cI^,^ 



0. 



where |0c) is the vacuum state with respect to operators c„^a- In the holomor- 
phic representation operators and Cn,\ are replaced by c-number functions 
^n,xiP) Zn,\iP) over which the integration is performed (see e.g. |^). One 
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can easily check it up by dividing the "time" interval [x^, x] into large num- 
ber of small ones and repeating step by step the derivation of the formula 
(55) inserting between the parentheses the unit operator 



I = I n<y^n,Ae-^-^<^^-^ \Z) {Z\ 



(51) 



where the state 



\Z) = exp lj2zn,XC^n,X \ \^c) 



has the properties 



Z) = Zn,x\Z) 

zS = H{zlz^)e-^2^^ 



{Z2 

H(c\c) 



and operator H is supposed to be normal ordered. 

The functional integral over cr(/3) gives A^„ = 1 in (|50D , i.e. the nonvan- 
ishing contribution to the matrix element comes only from the states where 
there is only one particle in each chain site n. Acting on these states hamil- 
tonian H transforms them into themselves since it does not change the total 
number of particles in the site. Finally, correlator (^81) is expressed through 
the matrix element in the kinematic where ~ x, ?/j ~ by the formula 



Ml , Ml V, A* AT, A* AT 



(Xi, ...,Xn) 



-< Oclll' 



TO, /I, 



,F(ct,c)e" 



:ln- 



-H 



'n,ii[ 



. 10. > . 



(52) 



The powers of z*{xi)zi{xi) in pre-exponential factor F{z*, z) (^91) do not con- 
tribute to functional integral (|48| ) (one can treat them as a result of insertion 
expression (|5T| ) between unit operators) therefore 

F(c*, c) = N |n'[i7i_i,i^(xi_i - Xi) + Hi^i+ie{xi+i 



Xi)] 



(53) 



where 



Hik - 
H 



n,n+l ; 



(54) 



prime has the same meaning as in Eq. ( ^Uj ) and means the normal ordering. 

Introducing helicity basis = (1, ±i) the two-body operators Hi^ act 
on the chain wavefunction ip as 



{H,k^p)++ = (H^k^p)— = 0, 



(55) 
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where all indexes except those at i, k positions are omitted. Introducing for 
each chain site an effective spin-| variable Sz and treating the positive helicity 
as = I state and the negative one as = — | the term Hik takes a form 

Tj f Ci ck I C* Ck \ f Qi Qk ^ 

-njfc — + — I o^o^ — - 

where 5'+ are the spin operators. Thus the total chain hamiltonian ( ^) 
turns out to be the hamiltonian of Ising model for one-dimensional magnets 
(for closed chain). 

However there is an uncertainity in this way of proceeding, namely, a 
general normalization constant that usually remains to be unknown in the 



functional integral approach. As a consequence hamiltonian H in (|38| ) is 
determined up to constant Eq which leads to the shift of the spectrum. To 
find Eq one has to compare the result with a case where the answer is known 
since Eq depends only on the way the functional integral is built up rather 
than a particular structure of the interaction. To this purpose the simple 
kernel 

can be taken instead of (0). This kernel is diagonal in polarization indexes, 
therefore the solution is given by expression (^). From the other hand it 
can be obtained through all functional integral machinery described above 
that leads to the quantum chain with hamiltonian 

n 

Due to the restriction Nn = 1 

H = Y.1 = N 

n 

and the asymptotic behaviour of the matrix element is the same as correlator 
( ^2D has. Thus one can conclude that energy shift Eq = 0. 

5. The solution of the Ising model is well known. After substituting the 
wavefunction of the Schrodinger equation 

the asymptotic behaviour of matrix element (152|) is defined by the maximal 
energy e^^J^. Then 



iy(i?) ~exp{l.W In^} 



and 



AN) ^ 



p(0 ~ exp J4£maa;^ln . 

V Xb 
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In the Ising model the zero-energy state (vacuum) is the state where all 
spins have the same orientation (e.g. when all they downward directed). An 
excited state with m turned-off spins is treated as the state with m quasi- 
particles (m < N) which are characterized by momenta ki. The excitation 
energy is given then by the formula 



e. 



(TV) 



Y^icoski + 1), (56) 



1=1 



where ki (or, more exactly, e*'^*) values are determined (for the closed chain 
case) by the algebraic equation of high order (see e.g. 0|) so there is not 
an analityc expression for them. Nevertheless Eq.(0) shows that e^^^, < 
2{N — 1) (the state with m = N obviously has zero energy). For large 
one can expect the maximum energy for m ~ y which leads to e^J^ = sN, 
where s ~ 1. Neglecting as running it gives for the structure function 



xbG2n{xb,Q ) ~ yo^j exp y 2sN ^ InQ'^ / ij,^\nl/xB. 

Thus the contribution of a single 2A^-particle cylinder is smaller than N 
Pomerons one (when Nc ^ oo Pomeron can be treated as a two-particle 
cylinder). Qualitative agreement of this asymptotics with formula (^2]) shows 
the incorporation of gluon polarization does not change the behaviour of the 
structure function. This is a consequense of small number gluons interactions 
allowed by cylinder topology. The stronger iV-dependence is to be expected 
for corrections. Therefore the investigation of the high twist evolution 
equation for finite Nc is of great interest. 
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Fig.l. The effective amplitude. 
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Fig. 2. An example of t-line amplitude. 
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